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We study the low-temperature low-frequency conductivity a of an interacting one dimensional elec- 
tron system in the presence of a periodic potential. The conductivity is strongly influenced by 
conservation laws, which, we argue, need be violated by at least two non-commuting Umklapp pro- 
cesses to render a finite. The resulting dynamics of the slow modes is studied within a memory 
matrix approach, and we find exponential increase as the temperature is lowered, a ~ (An) 2 e T ° // ' JVT ' 
close to commensurate filling M/N , An = n 

73.50.Bk,72.10.Bg,71.10.Pm 



M/N < 1, and a ~ e ( T o/ T ) 2/3 elsewhere. 
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The finite-temperature conductivity of a clean one- 
dimensional wire Q is a fundamental and much stud- 
ied question. Clearly the "bulk" conductivity of a wire 
in the absence of a periodic potential is infinite even at 
finite temperatures T. In this case the conductance is 
independent of the length of the wire and is determined 
by the contacts only. Surprisingly, much less is known 
about the conductivity in the presence of Umklapp scat- 
tering induced by a periodic potential. There is not even 
an agreement whether it is finite or infinite at finite tem- 
peratures for generic systems 0-0] • We shall show that 
the correct answer emerges when all relevant (weakly vi- 
olated) conservation laws are taken into account. Those 
conservation laws are exact at the Fermi surface and are 
violated by Umklapp terms away from it. We shall study 
the associated slow modes by means of a memory matrix 
formalism able to keep track of their dynamics. It will 
allow us to calculate reliably the low temperature, low 
frequency conductivity. 

The topology of the Fermi surface of a Id metal de- 
termines its low-energy excitations. Two well defined 
Fermi-points exist at momenta k = ±fcp , allowing us to 
define left and right moving excitations, to be described 
by ^L/R,a=u- We shall include in the fields momentum 
modes extending to the edge of the Brillouin zone, usu- 
ally omitted in treatments that concentrate on physics 
very close to the Fermi-surface. 

The Hamiltonian, including high energy processes, is 



H — Hll + H; 



EttU 



(1) 



Hll is the well-known Luttinger liquid Hamiltonian cap- 
turing the low energy behavior (!]], 



Hll = v F 
1 

~ 2 



p(xf 



dx 
2^ 



K v {d x e v f 



vp is the Fermi velocity, g > measures the strength of 
interactions, p = pl + Pb. is the sum of the left and right 



moving electron densities. In the second line we wrote 
the bosonized [jj] version of the Hamiltonian. Here v a , v p 
are the spin and charge velocities, and the interactions 
determine the Luttinger parameters K v with v v K v — vp, 
Vp/K p = v F + g/ir, v a /K a = v F - 9/tt. 

The high energy processes are captured in the subse- 
quent terms which are formally irrelevant at low energies 
(we consider only systems away from a Mott transition, 
i.e. away from half filling). Some of them, however, de- 
termine the low-frequency behavior of the conductivity at 
any finite T, since they induce the decay of the conserved 
modes of Hll (they are "dangerously irrelevant"). We 
classify these irrelevant terms with the help of two oper- 
ators which will play the central role in our discussion. 
The first one is the translation operator Pt of the right- 
and left-moving fields, the second one, Jo = Nr — Nl, 
is the difference of the number of right- and left-moving 
electrons, and is up to vf, the charge current of Hll ■ 

P T = Y / fdx (*L(-*&0*i* + *LH^)*l.) (2) 
J =N R -N L =Y1 I dx (*L*»r - *L*L CT ) (3) 

Both Pt and Jo are conserved by Hll, their importance 
for transport properties is due to the fact that both stay 
approximately conserved in any one dimensional metal 
(away from half filling): processes which change Jo are 
forbidden close to the Fermi surface by momentum con- 
servation. The linear combination Pq = Pt + kpJo can 
be identified with the total momentum of the full Hamil- 
tonian H and is therefore also approximately conserved. 

We proceed to the classification of the formally irrele- 
vant terms in the Hamiltonian. This classification allows 
us to select all those terms (actually few in number) that 
determine the current dynamics. H{ rT includes all terms 
in H — Hll which commute with both Pt and Jo, such 
as corrections due to the finite band curvature, due to 
finite-range interactions and similar terms. We will not 
need their explicit form. 
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The Umklapp terms H% m (n,m = 0, 1, ...) convert n 
right-movers to left-movers (and vice versa) picking up 
lattice momentum m2ir/a — mG, and do not commute 
with either Pp or Jq. Leading terms are of the form, 



Hk 



u 

9o,n 



(pl + Pr) 2 + h.c. 



(4) 



Km « dZrn £ / e^^ R ^ LaP - a + h.C. (5) 
Km » 9lm J e iAk ^*tf m tf Rl * Ll * L 1 + h.C. (6) 

with momentum transfer Afc„ jm = n2kp — mG. A pro- 
cess transfcring n > 1 electrons with total spin n s /2 
pointing in the z-dircction can be neatly expressed as 



n,m To \7> 
{2lTa) n 



,iAk n , m x e n/2(n0 p +n s 0<O _|_ /j c (7) 



a being a cut-off, of the order of the lattice spac- 
ing. In fermionic variables the integrand takes the form 
U^o'^R-iix+ja^^x+ja^L^x+ja^L^x+ja) 
(for n s = and even n). 

Note, though, that any single term H% m conserves a 
linear combination of Jq and Pp, 



[K m ,Ak nm J + 2nP T ] =0. 



(8) 



Indeed, a term of the form (0) would appear in a con- 
tinuum model without Umklapp scattering, but with a 
Fermi momentum hp = Ak nrn /(2n). In such a model, 
Ak nm Jo/ (2n) + Pt is the total momentum of the system 
and therefore conserved. The importance of this simple 
but essential conservation law has to our knowledge not 
been sufficiently realized in previous calculations of the 
conductivity. Due to this conservation law a single Umk- 
lapp term can never induce a finite conductivity! At least 
two independent Umklapp terms are required to lead to 
a complete decay of the current. Further, two incommen- 
surate Umklapp terms suffice to generate the rest. 

To calculate the conductivity it is necessary to keep 
track of the nearly conserved quantities and their relation 
to the current. We will develop a description of the slow- 
est variables using the Mori-Zwanzig memory functional 
j^,|]J^]. Approximations within this scheme amount to 
short-time expansions. In general, the short time decay 
of a quantity carries little information on its long-time 
behavior; this, however, is not the case for the slowest 
variables in the system, where the short time and hydro- 
dynamic behavior coincide. 

To set up the formalism j?j we define a scalar product 
(A|_B) in the space of operators, 



(9) 



(A(t)\B) = -J dX(A( t yB(iX)) 



where we use the usual Heisenberg picture with A(t) 

e 



which includes j\ — J, the full current operator. Stan- 
dard arguments lead to the electric conductivity, 



a(w,T) 



(10) 
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Here \ pq = (3{j p \j q ) is the matrix of the static j p j q sus- 
ceptibilities (as usually defined), and M is the matrix 
of memory functions given by the projected correlation 
functions of time-derivatives of the "slow" operators, 



HJq 



Q 



uj — QLQ 



dtjr (r L )rp- (11) 



The Liouville "super"-operator, L, is defined by LA = 
[H, A] and Q is the projection operator on the space per- 
pendicular to the slowly varying variables j p , 



(12) 



We assumed for simplicity that all j p have the same sig- 
nature under time reversal. 

The perturbative expansion of the memory matrix M 
is accompanied by factors 1/lu guaranteeing it is always 
valid at short times. It is also valid for small frequencies 
provided the slowly evolving degrees of freedom are pro- 
jected out (by the operator Q). Unlike the conductivity 
it is expected to be a smooth function of the coupling 
constants which can be perturbatively evaluated. 

We first consider a situation where some linear combi- 
nations of the j p are conserved by H, in which case an 
infinite conductivity is expected. We introduce V c , the 
projection operator on the space of conserved currents, 
and carry out the required matrix inversion to find, 



— > 0, T > 0) = a Icg (uj, T) + i 



(xic 1 x)n 

u> + iO 



(13) 



Ht Ae lHt . We choose a set "slow" operators ji,j2, ---jN 



where Xc 1 — 'Pc(PcXPc)~ 1 Pc- Within any simple (short- 
time) approximation, <7 rog (w,T) as defined above, is reg- 
ular (this approximation fails e.g. if some conserved cur- 
rent j is not included in Ji...jjv)- Hence the Drude weight 
D(T) is finite at finite temperatures, Recr(w — > 0) = 
2ttD(T)S(uj) = 7r(xxJ 1 x)ii^( w )- H is determined by the 
"overlap" of the physical current operator J with the 
conserved quantities %i s , s labeling the conserved cur- 
rents. Remarkably, our perturbative approximation is in 
accord with an exact inequality [^| for the Drude weight, 
D(T) > ^XXc 1 ^!!- Note that \ can b e calculated to 
an arbitrary degree of precision around a Luttinger liquid 
and that the lower bound can be improved by including 
more conserved quantities ^ . 

Now consider the more realistic situation where the 
previously conserved currents decay slowly (via Umk- 
lapp processes), in which case a finite conductivity is 
expected. We restrict ourselves to the two-dimensional 
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space spanned by v F Jo and P F , which we argue have the 
longest decay rate and dominate the transport. Here we 
approximate J ~ v f Jq to keep the presentation simple. 
This affects only the high frequency behavior of the con- 
ductivity There is a large number of other nearly 
conserved quantities. For example Hll + H^, the rele- 
vant low-energy model close to half filling, is integrable 
and therefore is characterized by an infinite number of 
conservation laws. We can, however, neglect them at low 
T if our initial model is not integrable, expecting that 
practically all conservation laws are destroyed by (for- 
mally irrelevant) terms close to the Fermi surface lead- 
ing to decay rates proportional to some power of T. This 
is to be compared to Jo and Pt which commute with 
all scattering processes at the Fermi surface, leading to 
exponentially large lifetimes. 

We now proceed to calculate the Memory matrix. To 
leading order in the perturbations we can replace L in 
(ID) by L LL = [H LL ,.] ||, since d t v F J and d t Pr are 
already linear in <?„ m . As Lll-Pt = LllJo — 0, there 
is no contribution from the projection operator Q. The 
memory matrix takes the form, 



M 



nm 



v 2 F (2n) 2 
— 2nv F Ak n 



-2nv F Ak nm 



where, 



/ 2v f /tt 
{ 



(4 



( 5 1) 2 M„(AU, W ) 



(14) 



Here F = [Jq, H^ m ]/(2n) (for simplicity we drop the in- 
dices n, to on F), and (F; F}° is the retarded correlation 
function of F calculated with respect to H FF . 

The memory function of the 4k F — Q process B^i 
was calculated by Giamarchi (|] , (not considering the ma- 
trix structure of M required by the conservation laws.) 
Higher Umklapps are considered in For n s = and 
even n the memory function due to the term (Q) can be 
analytically calculated, 



M n (Ak,u) 



2suv2nK r n 



2vraT 



4if"-2 



iS + , 1 - 2K n p )B(K n p - iS+, 1 - 2KI 



-B(K p - iS%, 1 - 2K n p )B(K n p - i$° , 1 - 2K n p )\ 



,,2-2n 



tt 2 T 2 (2K 1 })v p T 



aAk 



4K™-2 



,-v p Ak/(2T) 



where K n p = (n/2) 2 K p: B(x,y) = T(x)T(y)/T(x + y) and 



S± = (w±u p Afc)/(47rT), = S±{w = 0). The last line 
is valid for u) — and T <C v p Ak. 

The origin of the exponential factor is as follows: pro- 
cesses involving momentum transfer Ak are associated 



with initial and final states of energies u|Afc|/2, which are 
exponentially suppressed. If only charge degrees of free- 
dom are involved v = v p , otherwise v — min(w CT , v p ) = v a 
0. For T <c v a Ak nm , n s > and u> = 0, we have, 



, (15) 



M (Ak) ~ 

while for T » v p Ak nm : M n ~ T n2K " +n i K °~ 3 . 

Using the above expressions with only one Umklapp 
term leads to a finite Drude weight (cf eq([jT 



D(T) a - pKp 



1 



1 3(^Afc„ m V [ L + 4) 



1 



(16) 



in accord with the observation that one process H% m is 
not sufficient to degrade the current. 
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FIG. 1. The low frequency behavior of (t(u>) in the pres- 
ence of two Umklapp terms for two different T. The dashed 
lines are the result one obtains in conventional perturbation 
theory neglecting Q the matrix structure of M and the re- 
lated conservation laws. (#20 = 321 = 1, K p = 0.7, K a — 1.3, 
Ak 2 i = -1.5A/S20, thick lines T = 0.2, thin lines T = 0.18, 
ui and T measured in units of v p Ak2o). Note that two time 
scales appear - each describing the scale on which the asso- 
ciated conservation law is violated. The inset displays the T 
dependence of a(u> — 0). 

Only in the presence of a second incommensurate pro- 
cess HfL 1 is the dc conductivity finite, 



cr(T,W = 0) 



{Ak nm ) 2 /M n , m , + {Ak n , m ,) 2 /M n 



Tr 2 (nAk r , 



7i'A/c nm ) z 



(17) 



Note that the slowest process determines the low-T con- 
ductivity. The frequency and temperature dependence of 
the conductivity in the case of two competing Umklapp 
terms is shown in Fig. |l|. 

The commensurate situation Ak nm = requires extra 
considerations. Whether the dominant scattering pro- 
cess H nm will completely relax the current J depends 
according to ( [l3| ) on the overlap xjp t {[PT,H nm ] = 0). 
Using the continuity equation for the charge, xjp t can 
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be related to the deviation Ap = 2 An/a of the electron 
density from commensurate filling with the remarkable 
identity xjp t = 2An/a + 0(e- peF ). In a 3d lattice of Id 
wires, An is fixed by charge neutrality and is T indepen- 
dent, in a single wire with contacts An varies at low T 
with An(T) ~ T 2 /(mv 3 ), where the mass m is a measure 
of the breaking of particle-hole symmetry, e.g. due to a 
band-curvature k 2 /2m. In this case it is important to 
replace Ak nm = in Eqn. @ or @ by GAn(T). 




1/10 1/6 1/5 1/4 n 1/3 2/5 



FIG. 2. Schematic plot of log[cr] as a function 
of the filling n = 2/cf/G for various temperatures 
(vGf3 = 30,100,300,500,800) based on an asymptotic 
(eqn.(|l5|)) evaluation of (0). Near commensurate fillings n c 
the conductivity is strongly enhanced at low temperatures but 
drops at ti = n c . The inset displays the T-dependence of a 
for n — 1/3 and a filling very close to 1/3 (dashed line). 



typical ~ exp[c(/3«G) 2 / 3 ] (19) 

where c is a number depending logarithmically on T . At 
present we cannot rule out that various logarithmic cor- 
rections sum up to modify the power law in the exponent. 
We argue, however, that due to the exponential increase 
(U|) of a at commensurate fillings with exponents pro- 
portional to 1/Nq, the conductivity at small T at any 
incommensurate point is smaller than any exponential 
(but is larger than any power since any single process is 
exponentially suppressed). In Fig. || we show schemat- 
ically the conductivity as a function of filling becoming 
more and more "fractal-like " for lower T. 

Can the effects we predict be seen experimentally? The 
complicated structures as a function of filling shown in 
Fig. ^| are not observable in practice as they occur only at 
exponentially large conductivities. The T-dependence of 
the conductivity at intermediate temperatures, however, 
should be accessible, e.g. by comparing the conductivi- 
ties of clean wires of different length. Perhaps more im- 
portantly, it is straightforward to apply our method to a 
large number of other relevant situation, e.g. close to a 
Mott transition or in the presence of 3d phonons, as we 
will discuss in a forthcoming paper. 
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boldt Foundation and NSF grant DMR9632294 (A.R.). 



Which of the various scattering processes will even- 
tually dominate at lowest T? At intermediate tempera- 
tures, certainly low-order (small n) scattering events win, 
being less suppressed by Pauli blocking. At lower tem- 
perature the exponential factors in ( |j~5| ) prevail and the 
processes with the smallest Ak nm are favored. We first 
analyze the situation close to a commensurate point kp ks 
GM /(2Nq). The two dominant processes are H^ oMg 
with Ak NoMo w and H% lMi with Ak NlMl = ±G/N 
(or N X M = ±1 mod N ). The integer Ni of order N , 
N% = 71-/V0, depends strongly on the precise values of No 
and Mo. We thus find that the d.c. conductivity at low 
T is largest close to commensurate points with, 

a(k F » GM /(2N )) ~ (An(T)) 2 exp[/3«G/(2iV )] (18) 

but a ~ T~ N o K »- (No mod 2 ) 2i ^+ 3 if the density is exactly 
commensurate with \ An(T)\ < e -0Gv/(4N o ) _ 

To estimate the conductivity at a typical "incommen- 
surate" point or at commensurate points at temperatures 
not too low, we have to balance algebraic and exponen- 
tial suppression in (|l5| ) by minimizing —(3vG/(2N) + 
(7i7V) 2 if log[T] in a saddle-point approximation to the 
sum over all Umklapp processes in M, Up to logarithmic 
corrections we obtain iV^ ax ~ f3vG/(ji) 2 and therefore 
for a "typical" incommensurate filling, 
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